ITP-UU-10/39 
DFTT 27/2010 

D3-D7 Quark-Gluon Plasmas at Finite 

Baryon Density 

^ , Francesco Bigazzi Aldo L. Cotrone ^, Javier Mas 

^ ; Daniel Mayerson '^^ Javier Tarrfo ^. 

X) 

I ! " Dipartimento di Fisica e Astronomia, Universitd di Firenze and INFN Sezione di 

ly-^ ! Firenze; Via G. Sansone 1, 1-50019 Sesto Fiorentino (Firenze), Italy. 

* Dipartimento di Fisica, Universitd di Torino and INFN Sezione di Torino; 
Via P. Giuria 1, 1-10125 Torino, Italy. 
^ Departamento de Fisica de Particulas, Universidade de Santiago de Compostela and 
Oh| Instituto Galego de Fisica de Altas Enerxias (IGFAE); E-15782, Santiago de Gompostela, 

I Spain. 

Institute for theoretical physics, K. U. Leuven; Gelestijnenlaan 200D, B-3001 Leuven, 

■ Belgium. 

^ ■ ^ Institute for Theoretical Physics, Universiteit Utrecht, 3584 CE, Utrecht, The 

^ ■ Netherlands. 

■ bigazzi@fi.infn.it, cotrone@to.infn.it, javier.mas@usc.es, drm56@cam.ac.uk, 

l.j .tarriobarreiro@uu.nl 



X 



Abstract 

We present the string dual to SU{Nc) A/" = 4 SYM, coupled to Nf massless fundamen- 
tal flavors, at finite temperature and baryon density. The solution is determined by two 
dimensionless parameters, both depending on the 't Hooft coupling Xh at the scale set by 
the temperature T: eh ~ XhNf/Nc, weighting the backreaction of the flavor fields and 
6 ~ X'f^^^\b/{NfT^), where Hh is the baryon density. For small values of these two parame- 
ters the solution is given analytically up to second order. We study the thermodynamics of 
the system in the canonical and grand-canonical ensembles. We then analyze the energy loss 
of partons moving through the plasma, computing the jet quenching parameter and studying 
its dependence on the baryon density. Finally, we analyze certain "optical" properties of the 
plasma. The whole setup is generalized to non abelian strongly coupled plasmas engineered 
on D3-D7 systems with D3-branes placed at the tip of a generic singular Calabi-Yau cone. In 
all the cases, fundamental matter fields are introduced by means of homogeneously smeared 
D7-branes and the flavor symmetry group is thus a product of abelian factors. 
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1 Introduction and summary 



Heavy ion collision experiments at RHIC and LHC allow us to explore a relevant corner of 
the QCD phase diagram (high temperature and relatively small baryon chemical potential), 
where the theory is expected to be deconfined. Both the results collected during the ten- 
year run of RHIC [1] and the preliminary ones at LHC [2j actually indicate that a quark- 
gluon "fireball" is formed and behaves like a strongly coupled system: a liquid with very 
small viscosity over entropy density ratio. Holographic methods provide interesting tools 
to analyze these kind of systems. The simplest and best studied example is the conformal 
A/" = 4 SYM plasma which, unexpectedly, has proven to share some properties with the QCD 
one. This fact has stimulated further research works with the aim of refining this master 
holographic model, for example by adding fundamental matter fields. The latter has been 
performed mainly in the quenched approximation. 

In [3] some of the authors have presented a ten dimensional black-hole solution dual to 
the non conformal plasma of A/" = 4 SYM coupled to Nf ^ 1 massless fiavorsQ The 
latter were introduced by means of homogeneously smeared D7-branes O [6l [7], extended 
along the radial direction up to the black hole horizon. The smearing reduces the fiavor 
symmetry group to a product of abelian factors and allows a simple way to account for the 
backreaction of the D7-branes and thus to explore the "unquenched" regime in the dual field 



D3-D7 systems with D3-branes placed at the tip of a generic singular Calabi-Yau cone. In 
the zero temperature limit, the resulting backgrounds coincide with those found in [9j|f| 

In the present paper we extend the above construction to include a finite baryon density 
(or chemical potential, in the alternative thermodynamical ensemble) for the fiavor fields. 
Working on this problem with holographic techniques is especially interesting, taking also 
into account that there is no systematic way of dealing with finite baryon density in strongly 
coupled QCD (lattice QCD suffering from the so-called sign problem). 

We provide a novel gravity solution, dual to the above class of flavored plasmas in the 
planar limit at strong 't Hooft coupling. While the equations of motion we derive are 
completely general, the solution can be given in closed analytic form up to second order in 
eh ~ XhNf/Nc (where Xh is the 't Hooft coupling at the temperature T of the plasma) and 
6 ~ X;^^^\b/iNfT^), where rih is the baryon density. The gauge theories we focus on become 
pathological at some UV scale, developing a Landau pole. This is signaled, for example, by 
a running dilaton (accounting for the breaking of conformal invariance induced by the flavor 
flelds) blowing up at a flnite radial value. Correspondingly, the dual gravity solutions are 
not reliable close to that scale. Keeping eh small allows both to focus on a regime where the 

^AU the hydrodynamic transport coefficients of the model were derived in [4]. 
^For other holographic studies of thermal unquenched flavors, see [5]. 
Other solutions employing the smearing technique appear in [101 111) . 




The analysis was also generalized to A/" = 1 non abelian plasmas engineered on 
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solutions are reliable and to decouple the IR physics - which is the regime we focus on - from 
the pathological UV behavior. 

We also consider the regime of non-large baryon density, 5^1, both because it is the 
relevant regime for the RHIC and LHC experiments and because it allows to derive an 
analytic solution. Exploring the 5 ~ 1 regime requires a numerical analysis, that we plan to 
provide in the near future. 

The main results and the outline of the paper are as follows. In section [2] we present the 
action and the ansatz for the D3-D7 setup at finite baryon density. A set of second order 
differential equations is given in terms of the functions of the radial variable appearing in 
the ansatz. In section |3] we solve the equations analytically, in a perturbative expansion in 
e/i and 6 up to second order in both parameters. In section H] we perform the study of the 
thermodynamics of the system in the canonical and grand-canonical ensembles, checking the 
(non-trivial) closure of the various thermodynamic relations. We then explore the effects of 
the baryon number density on the energy loss of probes through the plasma, in particular 
on the jet quenching parameter. While the overall effect of flavors is to enhance the jet 
quenching [3], the effect of finite baryon density depends on the specific choice of comparison 
scheme of different theories. We finally provide some considerations on certain "optical" 
properties of the plasma, thinking about the possible gauging of the global U{1). Section [5] 
contains some concluding remarks. We also provide an appendix with some details on the 
ten-dimensional action, the equations of motion, the Bianchi identities and their solutions. 

The backgrounds we provide correspond to charged black holes in (slightly deformed) AdS, 
the charge being dual to a finite baryon density. The regime of validity is completely specified 
and the solution is totally reliable in that regime - there are no uncontrolled approximations. 
It is the first solution of this kind in the literature and thus it is suitable for the study of a 
number of physical effects of the baryon density. We hope to explore further the physics of 
this system in the future. 

2 Ansatz and effective Lagrangian 

The field theories we focus on are realized on the 4d intersection of Nc "color" D3 and Nj 
homogeneously smeared "flavor" D7-branes. The D3-branes are placed at the tip of a Calabi- 
Yau (CY) cone over a Sasaki-Einstein manifold X5, the latter being a U{1) fiber bundle over 
a four dimensional Kahler-Einstein (KE) base. The ambient spacetime, a product of 4d 
Minkowski and the CY cone, will be deformed by the backreaction of both kind of branes 
which respectively source a (self dual) F5 and a Fi RR field. As a result the lOd metric 
will be in the form of a warped product and there will be a running dilaton. Moreover, the 
backreaction of the D7-branes will induce a squashing between the KE base of the Sasaki- 
Einstein manifold and the fibration [9]. 
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Finite temperature is realized by placing a black hole in the center of the background 
[3]. The D7-branes extend along the radial direction up to the black hole horizon. Their 
embedding is described by a constant profile, implementing massless flavor fields in the 
dual gauge theories. In this work we are interested in switching on a chemical potential 
for the U{1)b baryon symmetry. The dual picture involves a non- vanishing profile for the 
temporal component At of the worldvolume gauge field on the D7-branes [12]. Through the 
Chern-Simons coupling, this field can source F3 and form fields. 

All in all we will be dealing with a general type IIB action given, in Einstein frame, by 



S 



where 



3 2 ^ 2^4^ 



+ S 



(2.1) 



Nf Nj 

— Jd7 ^ Jd7 



(2.2) 



is the contribution of the flavor D7-branes. The gravitational constant and D7-brane tension 
and charge are, in terms of string parameters 



T7 /i7 



1 



2kI, (27r)7^2c,/4 



(2.3) 



The smearing procedure [HI [6] amounts to a replacement 



D7 



(2.4) 



for any form Xg defined on the brane worldvolume Here ^2 is a form orthogonal to 
the individual location of the D7-branes. For an arbitrary Sasaki-Einstein space X5, it is 
proportional to the Kahler form J^e of the Kahler-Einstein 4d basis [H] 



gs^2 — —2QfJKE ■ 



(2.5) 



For massless flavors, is a constant encoding the density of D7-branes in the relative 
quotient space X5/X3 with X3 the subspace wrapped by each of the branes 



Qf 



Vol{X,)g,Nf 
4Vol(X5) 



(2.6) 



*The smearing of the DBI part of the flavor branes is described at length in [5]. 
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The equations of motion and Bianchi identities that follow from the action (12. ip . (12. 2p are 
given in appendix Rl 

Concerning the metric, we will consider the following ansatz, which includes a family of 
generalized squashed Sasaki-Einstein manifolds 

dslo = h-^/^[-bdt^ + dx'dx^] + h^'^[hS^F^da^ + S'^ds]^E + F^{dT + Ake?] , (2.7) 

where the Kahler two-form of the four dimensional base is given in terms of the connection 
one-form as Jke = dAxE/'^- The ansatz (12.71) contains two squashing functions F{a) and 
S{a) (with dimension of length), whose quotient F/S parameterizes the effect of the flavor 
backreaction. The dimensionless functions h{a) and h{a) account for the warping and the 
blackening of the spacetime, respectively. Thus, in particular, an ansatz with 6 = 1 is 
appropriate for the zero-temperature, uncharged solution. We have used the invariance under 
diffeomorphisms to choose a convenient holographic radial direction, a (with dimension of 
length""^), and as we will see, a — )■ — oo (0) in the IR (UV). Given the smearing procedure 
(12. 4p all the functions in our ansatz depend only on the radial variable. 

The finite baryon density is dual to a nontrivial worldvolume U{1) gauge field, 

= 2Tia A'^{(J) da ^ dt . (2.8) 

A consistent ansatz for the other fields is 

$ = $((7), ^2=0, F^=Q^{dT + AKE), F5 = Qe(l + *)V(X5) , (2.9) 

F3 = Fiasrfx^ A dx^ A dx^ - — — rft A fia + SQfJe'^^bF^ dt A da A {dr + Ake) ■ (2.10) 

S F 

In these expressions, V{X^) is the volume form of X^, is proportional to the number of 
colors 

Fi23 is a constant (of dimension length"^) which we will show to be related to the baryon 
density, whereas J = J{a) is a function (of dimension length^) that describes the effects of 
the backreactionj§ its contribution is dictated by the Cq potential 

Cq = J{a) dx^ A dx^ A dx^ A {dr + Ake) A ^2 , (2.12) 

which is the natural D5 charge sourced on the world-volume of the D7-branes by the gauge 
field through the last term in (l2.2p FI In [13] the system without backreacting flavors was 



^Obviously, J{(j) ^ Jke- 
"Notice that, with the smearing, 

^r J Mfi J Mth 



N, 
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studied, and an ansatz for was used that only contained the piece proportional to -Fi23- 
However, for our equations of motion to be consistent, we need the presence of the other 
components; thus, we see that J(cr) naturally contains the effects of the backreaction of the 
flavors. 

Inserting the whole ansatz into the lOd equations of motion and Bianchi identities one 
finally arrives at a system of equations which the reader can find in formulas (1A.13I) -( IA.19I) . 
It is possible to describe the whole system in terms of an effective one- dimensional action 
from which the equations of motion can be derived 

S^y^^[L,,,a, (2.13) 

where Vi^3 denotes the (infinite) integral over the Minkowski coordinates, and 
L,D = -^{hg'hf + 12{\og'SY + Slog' F log's -h'^ 

' ^ ^ log' h + 8 log' S + 2 log' F) - AQ}^ (2.14) 
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AbF^S' + 246^2^6 - -F^^.e^'bh^F'S' - -QybS' 



-4e'^/^FQfS^^-{2na'A',y + eH^F^S^ - 32Q}e~^bF^J^ - 8Q}{27ra' A',)J . 
The constraint equation (]A.20p is the zero energy condition H = for the Hamiltonian 

H = -L,D + Y.i^['^^ ^, = {6,/i,F,^,<|.,A,,J}. (2.15) 

Since the gauge field At enters only through its derivative it leads to a "constant of motion" . 
In principle this is a new free parameter which is related to the charge density. However the 
equations of motion link this to the value of F123 in the ansatz for F3 in ( 12. 9 p as we now 
show. Let us fix this constant of motion as follows 

1£ ^ 2WQ,Fi23 . (2.16) 



dA't 



Solving for A'^ it gives 



(Qci^i23 + 8QU)bFS^ 
2na'A', = . ' . (2.17) 

'lQQ)F^S^ + e-*(Q,Fi23 + 8Q)Jf 



Exactly the same expression is obtained from the equation of motion for the form field 
(see eq. (lA.lSP ). Thus, by enforcing the integration constant as in f l2.16p for consistency, we 
are putting the system partially on shell. On the other hand, this obscures the analysis when 
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it comes to computing the thermodynamical potentials holographically, since it means that 
the canonical momentum conjugate to A[ was already present in the original Lagrangian. 
We will comment on this later on. 

It is natural to use equation fl2.17p to eliminate A'^ in favor of F123 and this is usually done 
in one of two ways: obtaining the equations of motion from fl2.14p and then imposing eq. 
f l2.17p . or else, performing a Legendre transformation to the Lagrangian 



6L 



ID 



no 



no 



6A{ 



A' 



(2.18) 



and then taking the Euler-Lagrange equations from the transformed action. Either way the 
equations of motion coincide and are given by flA.lSp -f TA.lQp in appendix A. 



3 The perturbative solution 

In the uncharged case A'^ = F123 = J = 0, the following exact solutions for the functions b 
and h are readily found [3]: 6 = e'^^o'^, h = Qd^ — ^'^^°'^) / i'^f'o) ^ where tq is an integration 
constant of dimension of length. The black hole horizon is at a —t- —00 and the extremal 
limit is reached sending ro — )■ 0. In terms of a more standard radial coordinate r, defined 
in such a way that h = R^/r^ with = Qc/4:, one gets 6 = 1 — (rg/r)^ precisely as for 
the unflavored AdS^ black hole. The horizon radius Vh = tq is related to the temperature 
of the black hole. The whole solution in [3] also depends on the dimensionless combination 
e = Q/e*, which weighs the backreaction of the D7-branes and, in fact, can be read as a 
flavor-loop counting parameter in the dual fleld theory. 

Now, e runs as the dilaton and thus (as a common feature of backreacted D3-D7 setups) it 
blows up at a finite scale r^p (corresponding to a UV Landau pole in the dual field theory), 
rendering the supergravity approximation not reliable. Keeping e small requires restricting 
the validity range of our solution up to an arbitrary cutoff r* <^ r^p, such that 

VoKXa) . Nj 

e* = QfC * = -, — tA*^^ ^ 1 , 3.1) 

IGvrVoKXs) ' ^ ' 

where A* = Airgsc'^'Nc > 1, Nc,Nf > 1, $^ = $(rj and we have used 02. lip . In the 
uncharged case |3], this condition was used to find an analytic perturbative solution, up to 
order el, for the remaining functions S, F, $ appearing in the ansatz. The related integration 
constants were fixed requiring regularity at the horizon and matching with the T = solu- 
tion p] at the UV cutoff r^,. The resulting functions thus contained the dimensional (resp. 
dimensionless) parameters r/j,r* (resp. e*). The UV cutoff dependent terms resulted to be 
of the form of both power-like and logarithmic corrections. Formally sending the arbitrary 
cutoff scale r* to infinity, the first kind of corrections drops out, while the second kind can 
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be handled taking into account that the function e has a logarithmic running (accounting for 
the breaking of conformal invariance induced at the quantum level by the massless flavors) 
such that 



et = Q/e 



e* 1 + lot 



rh 



0{e 



(3.2) 



This procedure allows to decouple the IR physics from the UV one and to write down a set 
of solutions containing just rh and en as parameters. 

In the present charged case, we are going to follow the very same procedure. Here we 
have a further parameter to deal with: we will call it 5, and we will show that it is related 
to the dimensionless combination of temperature and baryon chemical potential (or charge 
density, depending on the thermodynamical ensemble). We will then derive an analytical 
perturbative solution taking both e* and 5 to be much smaller than one, deforming the finite 
temperature flavor backreacted solution obtained in |3]. 

As a first step, let us introduce a dimensional parameter 5 and consider the following 
redefinitions 



123 



5 



Qc 



J{a) 



3/2 



(3.3) 



The reason behind this choice will be clear in a moment. Inserting these expressions into 
( ]A.14[) - fA.19[) and rewriting the dilaton as $(cr) = + 0(cr), with 0((T*) = 0, one readily 
arrives at the following system of equations 



(log 6)" 
(log/i)" 
(log 5)" 
(logF)" 

(0)" 



e-'t'J' 



Y r ^o* 



+ ei6' Z , 
J'2 



+ ey -Z. 

* 2 



/2 



-Ql^^ + 2 ^ + 32 bF^J^ + 4 e ^^^^ 

lO. p3 clO 1 

-2hF^S^ + 66^2^6 - e, ^<i>l'^''-±^ _ 16 e^^^ e"* foF^ P _ ^^^^ 

4hF^S^ - ]^ey't> hS"" - 52 ^ + 16 el5^ e'^ bF^j' - 2 e"-^ 
-eJ -Z, 

1,2 p3 clO 

ey<P bS' + 2 e, e^^/^ + 2 e, e'^^'FS'Y - 32 e'^^ bF'j^ 



Y 

7/2 1 

I ^2 r2 -7 

+ ej -Z 



(1 + 8eJ)bFS^ 



+ 8e-^bF^J. 



(3.4) 



with 
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^ _ eHh'^F^S' 
The constraint equation (]A.20p reads 



(3.7) 




= -^log'/ilog'6 + ^(log'/i)2-12(log'^)2-41og'&log'^ 
- log' h log' F - 8 log' F log' S + ^0'^ 

4 T/2 1 \ 

-32te-^FV^ + - \z \ . (3.8) 

The system (I3.4p - (l3.8p allows for a systematic expansion of all the functions in powers series 
of e* and 5^. This is essentially the main effect of the scaling relations (13. 3p . Once all 
the functions have been solved for, the worldvolume gauge field can be obtained from the 
following relation 

which is already first order in 5. From this, we also deduce (as previously announced) that 
J(cr) takes the effects of the flavor backreaction into account. 

In order to integrate the system fl3.4l) - fl3.8p it is easier to switch to a radial coordinate y = 
e^'^o" with ro an arbitrary parameter of dimension of length. The dimensionless parameter 5 
referred to as above is then defined as 



where the factor 4 is introduced in order to make it precisely 5 = d oi ref. [12] (see also the 
equation (14. 7p in this paper). We keep the greek symbol, however, in order to stress that 
our parameter is going to be perturbatively small. 

Following the same procedure as in the uncharged case, we are now able to provide analytic 
solutions to the equations above, in a perturbative expansion in and 5. We will skip here 
the intermediate step where the solutions contain cutoff-dependent terms (this can be found 
in appendix |A] up to order e* 5"^ in eqs. (1A.22P to (1A.26P ) and focus on their effective IR 
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expressions. Introducing the IR parameter eo = ^('"o) ^-nd the radial variable r (defined again 
in such a way that the warp factor keeps the standard AdS form) , they read 



h{r) 
b{r) 



S{r) 



F{r) 



$(r) 



J(r) 
Mr) 



1 - - 

12 

1 + 

" 40 

320 
1 - 



2 - 







24 



5r; 

1 



17 

y 

^^0 



-log 



2 - 



1 + 



(3.11) 
(3.12) 



log(l 



1152 



log 

24 ^ r 



(3.13) 



3 1 



log 



1 + 



-33 + 22^ + 33 1 - 2^ log 



24+'° 



'0 

17 ^ 1 , ro 

log — 

1152 24 ^ r 



1 + ^ 



(3.14) 
(3.15) 



6o5^ 
' 40 

192 



^2 ^2 

22^ + 5-1 



3 1 



log 



1 + 



+ 2G'(r; 



-21 + 154^ - 35-^ + 21 1 - 2^ log 



$0 + eo log 

ro 48 



120 



26 



1-- 

^2 







1 + 3 log — ) log — - 3 Li2 
ro / r 



27r - 15^ + ( 11 + 18^ ) log 



4" 

^4 



14G(r) +... 



(3.16) 



14 log 2 



+G(r) + ...) , 



Ana' V 6 / V ^ / 



+ 



(3.17) 
(3.18) 



where G{r) — 27r^ 2-^1 |, 1, 1 — ^ j is an hypergeometric function and Li2{u) = Yl'^=i ^ 
is a polylogarithmic function. Notice that J enters at order Cq in the equations, hence only 
the leading contribution in J is relevant in the solution. 

The above solution must be supplemented with a Jacobian factor for the change of radial 



10 



coordinate Y{r) = da/dr which will show up in the coefficient of dr"^ 

eo (^rl{2r^ - rl){r^ + r^) - 2r<' log (^1 + ^ 



Y(r) = , } . + 



+ ^ (^-7rl + IQry - Mrlr^ + Mr"^ log (^1 + ^ 



+ ... (3.19) 



4 Physical properties of the dual plasmas 
4.1 Thermodynamics 

In the previous section we have derived a solution that, in essence, is a black hole dressed by a 
set of scalar (dilaton), vector (Maxwell), as well as higher rank tensor fields. Now we are going 
to extract the thermodynamical properties of the solution, providing, in turn, a ffist non 
trivial validity check of the latter by verifying the closure of the standard thermodynamical 
formulae. As in j3], all quantities are obtained in power series of our perturbative expansion 
parameters and, therefore, the relevant thermodynamic relations can only be verified up to 
the relevant order. 

To begin with, let us stress that tq does not coincide with the horizon radius rh- This 
radius is defined by h{rh) = + O(eQ, 5"^) and is perturbatively shifted by the baryon densitjill 

r, = ro (^1 + ^(1 - log2) - |^(11 - 171og2) + ... j . (4.1) 

Notice that > tq. Notice moreover that 

en ^ e(r,) = eo + 0{el) , e"- ^ e"^''-^ = e*" (^1 + - log 2)^ , (4.2) 

so that we can trade eo for eh, ^e*" for 5e*'", and so on, in all of our expressions. In particular, 
looking at f l3.18p . we get, to leading order in our expansion 



(4.3) 



which vanishes at the horizon as required to ensure IR regularity. 



^In fact, there are two radii r± that solve b{r±) = + 0{eQ, 6^). This is reminiscent of the case in the 
Reissner- Nordstrom black hole. The value presented here corresponds to the external radius = r+, i.e. 
the event horizon. Numerically we have checked that as S increases these two radii approach each other. 
However, a potentially extremal black hole r_ = r+ cannot be obtained within the range of validity of our 
solution, which is perturbative in S. 
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The temperature can be computed in the usual way giving 



T 



13 



384 



2 ~. 



13 



(4.4) 



The entropy density is derived from the horizon's 



TT 



2 Vol(X5 



7 

24 



(4.5) 



As for other thermodynamical variables, this expression, at first order in eh and 5^, precisely 
reproduces (for the case of massless flavors) the one found in the probe approximation in 
[12] . The 0{e1) terms are instead completely new. 

Concerning the charge density, a proper value is given by the integration constant in f l2.16p . 
as its definition coincides precisely with the electric field displacement. In terms of scaling 
constants we have, to leading order 



dA[ (27r)W 



(4.6) 



and making use of (12.61) as well as (14.41) this may be casted in the form 



TC 



7/2 _ 



Vol(X5)V2 ^ 

Using (13. ip with — )■ Vh this can be written as 



^hd 1 + -th 



(4.7) 



16 ^ ^ 



-Vol(X5)3/2 



(4. 



At leading order, and for the case Vol(S'^) = 277^ and Vo^S*^) = vr^, we exactly recover 
formula (A. 11) in |12]|§ In the expressions above, Uq is the quark density of the system, 
related to the baryon density by the number of colors riq = NcUb- 



4.2 Thermodynamical Potentials 

We proceed now to the calculation of the Helmholtz and Gibbs free energies, J-" and Q respec- 
tively!^ These can be either directly evaluated starting from the expression for the entropy 
density and using the standard thermodynamical relations, or they can be deduced holo- 
graphically. In the latter case they are identified with the (renormalized) on-shell boundary 

^Beware that \ours = 'ingsc'^Nc = 2\theirs- 

^The difference of these quantities with respect to the forms defined in section [2] should be clear from the 
context. 
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action for the gravity background, evaluated in the corresponding ensemble. Consistency of 
the solution requires that whatever method is chosen the results are the same0 

A quick look at the ID effective Lagrangian given in ( 12.14p reveals that on one hand At is 
a cyclic coordinate and, on the other hand, F123 is a Lagrange multiplier. Amusingly enough, 
the fact that the equations of motion for and At are consistent with one another imposes 
(I2.16p . which is nothing but the statement that, up to constant factors, F123 and At are 
canonically conjugate variables. Hence, as it stands. Lid contains both A[ and -F123, hence 
velocities and momenta. As a consequence the associated action corresponds to neither the 
canonical nor the grand-canonical ensemble. 



4.2.1 Canonical ensemble 

The Legendre transformed Lagrangian Lio given in (12.181) is the natural one to describe the 
system in the canonical ensemble, since it is fully expressed in terms of the baryon density 
parameter -F123. Therefore, evaluating the associated action on-shell we should obtain the 
Helmholtz free energy, J-'. As it is well-known, the action has to be supplemented with the 
standard Gibbons- Hawking term to deal with a well-posed variational problem. Even with 
this addition, the evaluation presents divergences which we deal with by subtracting the 
same quantity evaluated on the Euclidean solution at the same temperature but without a 
horizon and also with no chemical potential. This procedure is explained in appendix B of 
[3], where we refer the reader for details. We obtain the Helmholtz free energy density 



5 



V3 SVoKXs) ' 



which, consistently, satisfies the thermodynamic relation —df/dT = s. To check this relation 
it is very important to consider the dependence of and 6 on T. 

The logarithmic running of e = Q/e**^''^ (see also eq. (13. 2p ) and the map between the 
horizon radius and T give [3] 

^='f + 0{el~5'). (4.10) 
Moreover, for Xh{T) = 47r^,e*'^(^)A'c one gets by the same token 

^ = eH^ + 0{el~6^). (4.11) 

In the canonical ensemble we must keep the physical (dimensional) charge density invariant. 
From equation (14. 8 p we see that the dependence of 5 on T at fixed rib comes from solving as 

hT) = ^^,{l + 0{e,) + ...) , (4.12) 



^'^In the following we will not report the details of the holographic calculations, but we will just give a sketch 
of the needed ingredients. Needless to say, we have verified that the results deduced from thermodynamical 
relations agree with those found from holography. 
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with a a T- independent constant. Using (14.111) we obtain 

'd~5{T)\ 5 



dT 



(4.13) 



rib 



Using (I4.10p and (I4.13P it is straightforward to check that —df/dT = s with / and s given 
in (14.91) and (14. 5 p respectively. This is a strong proof of consistency. 

Next we can evaluate the ADM energy density of the plasma, again subtracting the con- 
tribution from the zero temperature and baryon density setup to get rid of the divergences. 
Following appendix B in [3] the final result is 

r5 _ r 1 / 1 . / 7~A1 

(4.14) 



E 



ADM 



8 Vol(X5 



l + ^e,(l + 252)+^e2 



4 



which satisfies the relation e = f + sT. From here, and taking again into account (I4.10p and 
(I4.13p . we obtain the heat capacity at fixed baryon number density 



de_\ 



TT 



2 Vol(X5 



(4.15) 



4.2.2 Grand-canonical ensemble 

We would like to obtain the thermodynamic quantities corresponding to the grand-canonical 
ensemble. It would be tempting to think that the correct Lagrangian density to use here is 
the original one Li£, in (I2.14p . However, this is not the case. As mentioned before, the fact 
that we have identified the canonical momentum conjugate to A'^ with F123 interferes with 
the Legendre transform, since we see that this parameter already appears in Lid. Were it 
not for this fact, the following inverse Legendre transform 



an 



Lid — 



dL 



ID 



123 



123 



Fl23{At) 



where Fi23{A^) comes from solving 



dL 



ID 



OF, 



(4.16) 



(4.17) 



123 



would bring us back to the original lagrangian Lid = Lid- However, notice the presence 
of the term with F^23 (I2.14p . which upon (14.160 will change sign. Therefore, the relevant 
Lagrangian for the computation of the Gibbs free energy is given in (I2.14p with a sign flip 
in the term with F^23- Doing that, and following the same steps as for /, we arrive at the 
result 



V3 



T7 = -P 



8 Vol(X5 



26' 



1 2 



(4.18) 
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where p is the pressure. Now we can compare this expression with the one for the energy 
density given in (14.141) . The interaction energy, given by 



c 



T4 16Vol(X5) 



-h > 



(4.19) 



is just the same as in the uncharged case. Hence, charge density of massless flavors does not 
contribute to the breaking of conformal invariance. This is not unexpected in field theory 
and, in a dual gravity picture, it is the same as what happens in the case of the Reissner- 
Nordstrom AdS black hole, where the presence of a net charge does not spoil the relation 
e = 3p. 

We now would like to check the thermodynamic relations in the present ensemble. Notice 
that the interpretation of 5 has changed and, henceforth, also its behavior with the tem- 
perature. The reason is that the physical parameter to fix in the grand-canonical ensemble 
is the dimensional chemical potential, /x. The difference among the two thermodynamical 
potentials gives 

vr^ ~ ( 1 



Using (14. 7p we can extract the chemical potential"] 

Now, the value of the world- volume gauge field in the UV (i.e. at r ^ r/j) is (see eq. (14. 3p ) 



Comparing the variation of the grand potential with respect to A 



t,uv 



.^V^^^J$^J^,j(,^l,A, (4.23) 



with the relation 

and taking into account formula (14. 7p . we get the connection between the UV value of the 
gauge field and the chemical potential 

At,uv =fi (l- ^eh] . (4.25) 



Using the particular value of Vol{X5) = tt^ we get the same expression at leading order in eh as in eq. 
(A.IO) in ^ after identifying p. = S/2. 
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Note that this differs at subleading order from the result Af^uv = /i as obtained in [12] ; this 
is thus an effect of the backreaction of the flavors. 

From ( 14.2ip we see that for fixed fi, 6 acquires now a dependence on T like 



~ 4 /vrVoKXs)/. / 5 



Using (14.101) again and (I4.26P we now get 

'd6\ 6 



.IT) r^l + f + O('0). (4.27) 



With this scahng (and equation (I4.10p ) it is easy to check that the thermodynamic relation 
—duj/dT = s holds at the required order. Using again (I4.27P we can obtain readily the heat 
capacity at fixed chemical potential 



(4.28) 



from which we can extract the speed of sound 

The parameter 6c = l — Sc^ is related to the breaking of conformality as before. As it happens 
for the interaction energy, it does not receive corrections from the presence of finite baryon 
density. Using the arguments in [3], this suggests that the bulk viscosity is not affected by 
the presence of a finite baryon density on the systemic 

Susceptibilities 

The "quark" susceptibility!^ is 



X 



^^^^ NjNc T^(l + -eh] . (4.30) 
4Vol(X5) ^ V 6 7 ^ ^ 



At leading order this matches with the result obtained in the probe approximation in [14j 
for the flavored A/" = 4 SYM case (where Vo^Xg) = 271^ , Vo^Xg) = tt^). The other three 



^^There is some mismatch between our results for the speed of sound and the heat capacities and cor- 
responding resuhs found in the hterature in the probe approximation. The precise closure of the thermo- 
dynamical relations and the consistent independence of the conformality breaking effects from the baryon 
chemical potential in the massless flavored case, let us be confident of the correctness of our results. 

^■^As usual, "quark" is an abuse of language for "fundamental matter" . 
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susceptibilities are 

ff^u ds Svr^ ^ „ „ / 1 11 2\ 7rVol(X3) o A 1 ^ 

(4.31) 

and 



iV/iVe /i T (^1 + ^e;,^ . (4.32) 



d^&T dTdfi 2Vol(X5 
The determinant of the susceptibihty matrix equals xC where 

The second term is parameterically smaller than the first one, thus the theory is thermody- 
namically stable. 



4.3 Probe parton energy loss 

In order to estimate how the finite charge density or chemical potential influences the energy 
loss of an energetic probe parton traveling through the plasma, we will make use of the 
results in [151 IIS] • f his approach, the parton looses energy through bremsstrahlung due 
to its interactions with the strongly coupled medium. The amount of energy loss, which is 
ultimately the cause of the jet quenching in the strongly coupled plasma, is encoded in a 
transport coefficient termed g, the "jet- quenching parameter". The latter can be derived, 
using the eikonal approximation at high energy, as the coefficient of in an almost light-like 
Wilson loop with dimensions L~ ^ L. The Wilson loop is easily calculated in string theory. 

In [17] a formula to extract the jet quenching parameter from a general dual gravity 
background was derived!^ On our solution, rewritten in the variable y = exp(4rQ(T), the 
formula gives 



q-'=na' fe-t ^=dy = n [' ^-^.uIlhFS'dy , (4.34) 

Jo gxxVQxx + gtt Jo Ar^ye- V 1 - 

where in particular g^x is the metric coefficient in the spatial direction involved in the tra- 
jectory of the parton. For our purposes, it will be enough to analyze the result up to order 
e/i5^, which reads 



""'^^^^^ - ^ ■ - ' ■ (4.35) 



VvoKx^r 



4^ rp3 



2 + TT / ~, 



* We took into account a different factor of \/2 between the definitions in [TBI and [17] • 
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with a positive constant c = 0.8675650 
4.3.1 Comparison schemes 

Formula f l4.35p . without any prescription to compare different theories, would imply that 
a finite charge density (or chemical potential) increases the jet quenching. This conclusion 
would depend on considering the jet quenching parameter in theories with different numbers 
of degrees of freedom. From a phenomenological perspective, it is more useful to compare 
theories keeping the number of degrees of freedom fixed. In order to fully appreciate the 
effects of flavors at finite charge density, we will compare our flavored theory with the un- 
flavored M = A SYM. 

For an estimate of the number of degrees of freedom, the most used observables are the 
entropy density (14.51) . and the energy density (14.141) . In large Nc theories, it makes sense to 
keep one of these two quantities fixed by either fixing the number of colors Nc and changing 
the temperature T, or the other way around. First, let us consider the fixed N^.-, varying T 
comparison scheme, used for example in [18]. For the sake of definitiveness, we also keep 
the energy density fixed; the qualitative result which will follow is unchanged if we fix the 
entropy density instead. We get 



1 _ ( 1 + 25^ 



(4.37) 



where the subindex refers to either the /favored or the unfiavored A/" = 4 theory. Plugging 
this result in formula (I4.35P we read 




TT — 1 / ~9 

-eh f 1 - 0.719 5^ 



(4.38) 



where both 't Hooft couplings Aj,Au are evaluated at the horizon. Since N^. is fixedH if 
we keep the coupling as fixed, we obtain the result that while the jet quenching parameter 
is enhanced by the addition of flavors [3], the finite charge density (or chemical potential) 
actually reduces the enhancement 



This qualitative result is unchanged if, following strictly [18] , we also allow for the variation 
of Us- The latter can be adjusted in such a way that the force between two external quarks 
at the screening length Lc,n of the unflavored plasma, aqq = 3L^„y(Lc,u)/4, is kept fixed. 
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The constant c is related to the resuh of the integral 



-(8 - Syi^ + 2y)(l -y) + x/T^(3 + 2y)(2arctanh[,/T^ + logy) ^ 
^^Remember that A ~ asNc- 

^^Remember that the term in ChO^ is by definition a perturbation of the one in e^. 
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The potential V{L) is calculated numerically by standard formulas, reported for example 
in [19]. We calculate aqq for different values of Nf (equivalently, et) and 6. At zero charge 
density, as we increase Nf, aqq increases, so in order to keep it fixed the coupling A/ must be 
But the effect is small, so as anticipated we get in fl4.38p that the jet quenching 
parameter is enhanced by the addition of flavors. Along the same lines, if we switch on the 
charge density, we find that this enhancement is reduced. 

We can now analyze what happens if we use the alternative comparison scheme introduced 
in [3], where the temperature T is kept fixed and the number of degrees of freedom is kept 
constant by varying the number of colors Nc- By fixing the energy density (I4.14P we get 



Thus, while again the overall effect of the flavors is to enhance the jet quenching, in this case 
the charge density actually increases the enhancement. The result is qualitatively the same 
if we keep fixed the entropy density (14. 5 p instead of the energy density. 

To summarize, while in all the cases the effect of fundamental flavors is to enhance the 
jet quenching [3], the net effect of a finite charge density (or chemical potential) actually 
depends on how we compare different theories. Namely, we found that the charge density 
reduces (resp. increases) the enhancement of the jet quenching due to flavors in the varying 
T, fixed Nc (resp. fixed T, varying A^^) comparison scheme. We currently have no intuitive 
explanation for this behavior. 

4.4 Remarks on "optical" properties 

The solution presented in this paper is dual, in the hydrodynamic limit, to a charged rela- 
tivistic fluid. Recently, it has been observed in |20l|2T] that, generically, this kind of fluids has 
interesting uncommon optical properties, such as negative refractive index, and exhibits ex- 
otic phenomena, like the propagation of additional light waves in certain frequency regimes. 
It is understood that one can speak of "propagation of light" by thinking of gauging the 
global U{1) charge: the gravity solution allows anyway for the calculation of the current 
correlators involved in the optical analysis. Alternatively, one could maybe think about 
"optical properties" of baryon charge waves (e.g. dispersion relations). 

^®This corrects a statement in footnote 20 of [3]. 




(4.39) 



Keeping fixed also the coupling constant for simplicity, we obtain 




(4.40) 
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To be concrete, in [20] it was shown that every relativistic fluid with finite charge density p 
{uq in our solution), having a dispersive pole in the (transverse) current retarded correlatocj 

Gt{uj, k) ~ + i^(0, k) , (4.41) 

displays negative refraction, i.e. the phase and group velocities have opposite directions, 
for small enough frequencies cu^ < 47rg^i3, where q is the U{1) gauge coupling, if the term 
P(0, k) can be discarded. The latter requirement was verified by explicit gravity calculations 
in [201 ini [22]. In formula fOTj) one has 

S = -^, ^ = ^' (4-42) 

e+p e+p 

where rj is the shear viscosity of the fiuid. In the case at hanc@ 

^ ^Vol(X3) , ~, V=^-'-^(l + . (4.43) 

Thus, assuming again that the term P(0, k) does not affect the result, the charged D3-D7 
Plasmas would exhibit negative refraction for 

Actually, in a medium with large spatial dispersion, there is the possibility of the propaga- 
tion of two "light waves" for each incident wave, due to the fact that the dispersion relation 
admits two different solutiona^ 



This effect is due to the coupling of the incident wave with "exciton" quasi-particles in 
the medium, which in the dual gravity description are quasi-normal modes. One of the 
two indexes should correspond to negative refraction for small u, while the other one is 
connected to the existence of an additional light wave. The propagation of both waves in 
certain frequency regimes could be inferred from the study of an "effective refractive index" 

Ueff = . (4.46) 

ni+n2 ^ ^ 



^^Clearly, here w, k are the frequency and wave vector. 
^°We make use of the relation rj/s — 1/4% [23] . 
^^Again, assuming that P(0, k) is ineffective. 
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Re[n] 




Figure 1: The real part of the two refractive indexes rii (upper hue) and n2 (bottom hne), 
with the effective one ne// (central line), as functions of the frequency (we are using unit 
temperature). The plot refers to the A/" = 4 theory with the following parameter choice: 
= 100, Nf = 20, A;, = 50, g = 1, 5 = 0.1. 

The results for the real part of the two refractive indexes and the effective one are summa- 
rized in figure fTll^ One of the two indexes displays negative refraction for co < coc, with ujc 
given precisely by f l4.44p F^ The interpolation of the effective index Ueff from one refractive 
index at small frequencies to the other one at large frequencies, which can be seen in the 
plot, signals that in the intermediate regime both light waves actually propagate. 

In conclusion, the outcome has very similar features to the ones found in the analysis of 
the RN black hole cases in [Sni EH 122]. Thus, we seem to see the expected features already 
at the level of hydrodynamics. It would be very interesting to study the Green functions 
directly in the gravity setting to understand whether this behavior is confirmed and to what 
extent the hydrodynamic approximation is reliable. 

5 Conclusions 

Following the strategy initiated in [3|, in the present paper we give the crank another turn, 
by introducing a new deformation that encodes holographically the backreaction of charged 

^^We thank A. Amariti, D. Forcella and A. Mariotti for their crucial observations on the results. 

^■^The imaginary part of one of the two indexes, Im[n2], would be negative for lo < ojc- This means that 
we have to choose, for oj < oJc, the negative branch of the square root of in (j4.45|) . which results in the 
negative refraction. We also find a particular value of uip > uJc for which the two indexes exchange (i.e. the 
continuous functions are made by gluing the two indexes aX Up). This is probably an artifact of Mathematica. 
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degrees of freedom at finite density. Following the same logic we provide an analytic solution 
which is perturbative in the deformation parameters, and explicitly write it down up to 
second order in both eh and 6. Moreover it is an effective solution. This means that only 
IR quantities can be reliably computed from it. In this case, the IR energy scale, /x/r, is 
set by the temperature, and all results have to be understood as being correct up to terms 
of order Hir/ Huv = T/ fiuv where fiuv is an ultraviolet cutoff energy scale above which the 
ultraviolet completion sets in. Still a lot of phenomenologically interesting prediction can be 
made both for thermodynamic quantities as well as for transport coefficients. We have shown 
that the thermodynamics is consistent, and we are in a stable phase. So far, nothing can 
be said about phase transitions. The solution contains two horizons, and has the potential 
to exhibit an extremal solution that would allow to explore the zero temperature - finite 
chemical potential axis. However such a systematic search of the phase space in the (T, /i) 
plane lies outside the domain of validity of the analytic solution presented here. It would 
require the use of numerical techniques, along the lines of the recent work [21] and certainly 
constitutes a natural continuation of the present work. 

Concerning transport phenomena, we have studied the jet quenching parameter and found 
that the effect of the net baryon density is positive or negative depending on the way one 
compares different theories with and without flavor. In the present context this effect is 
subdominant as compared to the one driven by the presence of flavor. Again, the fate of 
this at higher finite values of the charge density can only be established numerically, and its 
importance is evident from the fact that this is one of the few windows to phenomenology 
for a rather wide class of models. It is also worth comparing this result with the information 
one can obtain from the evaluation of the drag force. 

Last but not least, we have analyzed the possibility of exotic optical phenomena in the 
case that the global U{1) could be gauged and treated as the electromagnetic one. Following 
[201 l2Tj , we have checked the presence of a negative refractive index and the propagation of 
additional light waves in certain frequency regimes. Our analysis concerns the expectations 
for the optical properties, based on hydrodynamic considerations. It would be very interest- 
ing to go beyond the hydrodynamic approximation by the explicit gravity calculation of the 
relevant retarded Green functions. 
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A Technical details 

A.l Equations of motion from the ten dimensional action 

The equations of motion derived from the action S as given in (12. ip and (12. 2p are [H 



Rmn - ^9mnR = ^ (Om^On^ - ^gMNdp^d^<A + ^e^* (fimFin - ^gAmF^ 



^ T7 77 PQRS 2k, SSfi /\ \ 

+ ^F5MPQRSF^^ +-j=j^, (A.l) 

d(e2**Fi) = -e^H.A^F.-^J^'AQ^, (A.3) 

d{e'^*F,) = -H,AF, + ]:T^AQ2, (A.4) 

d{*F,) = dF, = HsAFs-^T^An2, (A.5) 



d{e~'^*H^) = e*Fi A *Fs - F5 A F3 + e*-^^- det(^ + e-^/^r)6^^\D7) , (A.6) 

Oj- V 

where the last term in the equation for has to be meant as an eight-form: in particular 
5^'^\D7) is a short-hand notation for the form which arises taking the derivative w.r.t. J-" of 



^■^In our conventions = ^{FpY^°-^---°-J'{Fp)a^a2...ap (also for H^). Notice moreover that self-duality of 
-Fs has been imposed. 
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the smeared DBI action, where a/— det(^ + e~'^/'^J-')d^x is replaced by 



det(^ + e-*/2J^)|fi2M^°a;. (A. 7) 

The expression for the modulus of the Q2 form will be given in a moment. 
The Bianchi identities read 

dF, = -gsn2, (A.8) 
dFs = /J3 AFi - J^A^]2, (A.9) 
dHs = 0. (A. 10) 

Finally, the Bianchi identity and EOM for the brane field J-" are 

dJ^ = H3, (A.ll) 

d (e*^y/- det(^ + e-^/^J^)6'^'\D7)^ = d{. ■ ■) , (A.12) 

where the dots represent the terms from the WZ part of the flavor action. 

After inserting the ansatz introduced in (12. 9p the equations of motion ( 1A.3I) (I A. 4 1) and ( lA.SI) 
as well as the Bianchi identities (]A.8|) and (]A.10|) are automatically satisfied. The l.h.s. of the 
equation of motion (1A.6P for vanishes identically upon imposing the ansatz, and from its 
right hand side (explicitly rewritten taking into account eq. (1A.7I1 . with |f22| = ^Qfh~^^'^S~'^) 
we can solve for A[ to obtain 

27ia'A[ = ' . (A.13) 

IQQ)F^S^ + e-*(QA3 + 8Q)J? 

This is an important relation, and interestingly enough, it automatically solves the equation 
of motion for At. This is because flA.12p is precisely the radial derivative of flA.6p . The 
Bianchi identity (lA.lip is identically null. 

The Bianchi identity for F3 f lA.9|) gives a nontrivial relation between A'^ and J . Using 
(IA.13P this can be expressed as an independent equation for J which we have added to the 
list below as eq. flA.lQp . 

Finally we come to the Einstein-dilaton equations of motion (lA.ip and flA.2p 

(log^r = 4Q;^ + 64Qje-nFV2 + 8Qje-*-^ + Z, (A.14) 
(log/.)" = -Ql^,+2Qf^ + ?>2Q)e-HF'j' + AQ)e'''^^ + \z, (A.15) 

7,2 p3 clO 1 

(log 5)" = -2hF'S' + QhF^S^ - Qi^-^l^ y - leg^e-'^ - , (A.16) 
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(logF)" 



1 Y ^2 



4fcF*>5^ - ^qy^ bS"" - Qf— + IGQje-'^bF'J' - 2Qje 



T/2 



(A.17) 



l2 ClO 1 



^4^2 



(QA3 + 8Q2J)6F54 



^i6Q}F^s^ + e-*(geFi23 + sgjj)^ 



2 

(A.18) 
(A.19) 



and the constraint that comes from fixing the reparameterization gauge invariance in the 
radial variable reads 



= 



log' h log' b + ^(log' hf - 12 (log' Sf-4: log' b log' S 
- log' 6 log' F-8 log' F log' S + ^$'2 



4^2g3*/2^3g^^l0 



*n2 7/2 



In these expressions we defined 
X = (27r«'>i;)2e*/2F^2 . y ^ y-(27rQ;'A02 + b'^e^F'^S^ ; Z = F^^sgC* fe/i'F^^^ . (A.21) 



(A.20) 



A. 2 Solution with explicit UV cutoff 

The solution with explicit dependence on the position of the radial cutoff r* to first order in 
e* reads 



b{r) 
S{r) 



^2\ ^2 



-log 



1 + ^ 



1 + ^ 



+ 0{el) , (A.22) 



= 1 + 



36 2rt - 



' 40 



f 6 - 44) + 4 (12r4 - 6r^) + (2rV2 - Gr^) + 4 (r^ - 2r^) 

\ ' / ' ' * 



'2rt - rl 



1-2^) log 



1 + 



1 + ^ 

' * 



+ 0{e 



(A.23) 
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F{r) ^ ^ 3rf + r4-2r^ 







36 2rt - 4 



2. /rf ( 6 - 444 + lO^i) + ( 124 - 6^') - f 6 - 22^ + 5^ 



40 V 2r4-r^ 

Tn 2r^ — rf. 16rf — Art Gn(r) ( , , 

^ * , " ^77^ - 3 1 - 2^ log 



1 + 



1 + 



+0{et) , (A.24) 



r 



$(r) = ^, + eJog{-]+0{ei), (A.25) 



S 4 4 



Jir) = -^ZTT^j + ^i^*)^ (A-26) 







with Go(r) = 27r^ 2^1 (|, |, 1, 1 - ^) . 
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